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^ • Abstract 

0> ■ Chiral symmetry at finite temperature is studied using the Schwinger-Dyson equa- 

tion. We calculate numerically the critical temperature using the Schwinger-Dyson 

lO ' equation with the gauge parameter that depends on an external momentum. The crit- 

O ' 

O ! ical temperature obtained by this method is similar to that with the Landau gauge 

and wave function renormalization constant 1. Moreover, the gauge invariance in the 

ladder approximation is examined using our method. 
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§1. Introduction 

In quantum chromodynaniics (QCD), chiral symmetry is approximately realized at the 
lagrangian level, and it is broken by the strong interaction at zero temperature. On the other 
hand, at finite temperature and density, the broken symmetries are restored. Therefore, one 
expects that finite temperature and density QCD have various phases. For example, the 
quark-gluon plasma (QGP), where quarks and gluons are deconfined has been observed at 
the Relativistic Heavy Ion Collider (RHIC).-'' The QGP is examined at the Large Hadron 
Collider (LHC).EJ 

To study the chiral symmetry breaking, we need a nonperturbative treatment due to the 
strong interaction. As a theoretical approach to understand the phase structure of QCD, 
the lattice QCD simulation is a powerful method.l^ The QCD phase structure at finite 
temperature is extensively studied by this method. The lattice QCD simulation derives the 
results that are consistent with the experiment at RHIC.® 

Another theoretical approach is the Schwinger-Dyson equation (SDE) method.^' The 
SDE is a valid method for both finite temperature and density. The lattice QCD for large 
chemical potentials is still inadequate. Moreover, the numerical analysis of the SDE does 
not require a large-scale computer like in the case of lattice simulation. When we solve the 
SDE, we need appropriate approximations, because the SDE is a group of infinitely coupled 
equations. 

Although the ladder approximation is usually used to solve the SDE, the SDE with this 
approximation depends on a gauge parameter. Owing to the gauge parameter dependence of 
the SDE, observable quantities, which should essentially be gauge-parameter- independent, 
depend on a gauge parameter. However, at zero temperature and density, the ladder ap- 
proximation Ward-Takahashi identity (WTI) is guaranteed by choosing the Landau gauge. 
Thus, the Landau gauge is used at zero temperature and density.'^ 

By contrast, at finite temperature and/or density, the ladder approximation WTI is not 
guaranteed trivially.'^ For this reason, there are no gauge parameters that have a clear 
advantage at finite temperature and/or density. The Landau gauge is adopted from the 
analogy of zero temperatureSI'''Si' . On the other hand, the Feynman gauge is adopted for 
convenience.'^ 

The purpose of this paper is to perform the numerical calculation of the ladder approxima- 
tion SDE to satisfy the WTI at finite temperature. For this purpose, we employ the method 
with a gauge parameter that depends on an external momentum. In this method, the WTI 
is satisfied using this functional gauge parameter. Such an idea, in which a gauge parameter 
is treated as a function, was studied at zero temperature^ and was used in QED with real 



time formalism at finite temperature.'^ (The formulation at finite temperature is shown, 
e.g., in Ref. 12).) However, the numerical calculation of the SDE in real time formalism 
uses further approximation (IE approximation®) in addition to the ladder approximation. 
Moreover, the results obtained with real time formalism using IE approximation at the zero- 
temperature limit do not correspond to zero temperature.'^ Thus, the present SDE in real 
time formalism is insufficient. Hence, in this paper, we use the method with the functional 
gauge parameter in QCD with imaginary time formalism. Since the numerical calculation 
of the SDE in imaginary time formalism need not use other approximations, this formalism 
is more reliable. 

The paper is organized as follows. In ^ we review the zero-temperature and finite- 
temperature SDE. In ^ we provide numerical results with the functional gauge parameter. 
The numerical method for solving the SDE is iteration. Then, we calculate the critical 
temperature. A summary and discussion are found in §4. 

§2. Schwinger-Dyson equation 

2.1. Zero-temperature SDE 

The SDE is derived using the CJT effective potential.'^ The SDE for quark is given by 

G-\p) = S-\p) - zg'C, I ^^^^D^-{p - q)G{q)r,{p, q), (2-1) 

where C2 is the Casimir operator, S{p) is the free quark propagator, D'^'^ is the exact gluon 
propagator, /^ is the exact quark-gluon-quark vertex, and G{p) is the exact quark propaga- 
tor. 

We use the chiral limit for the free quark propagator. In the ladder approximation, '^ the 
SDE is written as 

G-\p) = S-\p) - zg'G, I ^.l.D'o^iP ' Q)Giqhu, (2-3) 

where Dq"^ is the free gluon propagator, 
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Here, ^ is a gauge parameter. From this equation, the SDE is divided into coupled equations 
for scalar functions A{p) and B{p). 

In QCD, we use the improved ladder approximatiorP^ in which a coupling constant is 
replaced by the running coupling constant gip^iQ^)- The running coupling is included in a 
momentum integral. 

After performing angular integral,'^''^ 

^^"^ IQ-kH'^ J, A\s)s + B\s)^ IQtx^ J I A\s)s + B^s)' ^ ' 

^^WM^+O /■' sB{s) , ^2(3 + r. 9\s)B{s) 

^'~ IQ-kH J^ A\s)s + B\s)^ 167r2 / 'A\s)s + B\sy ^ ' 

where / = p\ and s = q"^, the index E shows four vectors in Euclidean space. In the 
ladder approximation, the wave function renormalization constant must be unity to satisfy 
the WTI, that is, A{1) = 1. This equation for A{1) shows that A{1) = 1 if the Landau gauge 
^ = is adopted. Thus, the WTI is satisfied using the Landau gauge at zero temperature. 

Although the WTI is extended to the Slavnov- Taylor identities in QCD, the Slavnov- 
Taylor identities become the WTI type to omit ghost-quark scattering kernel.'^ Since the 
ladder approximation corresponds to this situation, we use QCD with the WTI type. 

2.2. Finite-temperature SDE 

The Feynman rules in imaginary time formalism are summarized as follows: 

free quark propagator : S(p) = (po = iUn = 2niT(n -\ — 

7^p^ — m \ V 2 

tree giuon propagator : D^ [k) = — h q, ^ [ko = ^^i = ^^'^-^ ^j 

(27r)V(p+---) ^ -^(2vr)X...53(P+---), 

Here, a;„ is the Matsubara frequency. By this replacement, the improved ladder approxima- 
tion SDE for quark in imaginary time formalism is given by 

"^'^ -'' ^' "''- -"'{p-q)G{qh., (2-5) 






where p^ = (2i7rT(n + 1/2), pi) and q^ = (2i7rT(m + 1/2), qi). If we use the exact quark 
propagator form, Eq. fl2-2p . the coupled equations for A'^{x) and Bn{x) are written aCj 
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where, x = |p| and y = \q\. Alternatively, if we use the general form of the quark propagator 
at finite temperature, 

the coupled equations for A^(x), -B^(x), and C^(x) are 

C^n(x) = 1 + ^^ $: / dyg\-p\ -q')y 

-Cra{y){h + h) - AUy)h + i{Cra{y)h + Am{y)h) 



(2-8a) 



X 



Cl{y)ql - A1{y)y^ - BlXy) 



Mx) = l-^,J2 dyg\-p'^-^')y 

m " 

-C^(y)gi + AUy)iH, - Hs) + acUy)Hi + AUy)H,) 
"" cUyH - Aliy)y' - siiy) 

R r ^-^rr^llx- r ri 9\-p\-q')yBm{y) 
"^ ^ ^ 87r% 4^ io ^C7^(l/)gg - Al{y)y- - 5^(|/) 

, (Po - go)^ - (x + i/)2 
X log ■ 



(2-8b) 



(2-8c) 



(Po-go)^ - {x-y)'^' 

The explicit expressions of L, I, and if are given in Appendix A. Equations fl2-8al) - fl2-8cp 
are the general coupled equations derived from the ladder approximation SDE at finite 
temperature. Those four scalar functions have the relation for n, e.g., Bn{x) = -B_„_i(x). 



*) Since a propagator at finite temperature is Euclidean, Eq. (|2-2p is replaced by Gn{x) 
-1/(A;(x)p^7'^-B„(x)). 
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In this paper, we use the following form for the running coupling]^ 

9'{-P^ -1') = TT#L_ X ^ i + "-at^;;:S°'' . ' tj<t< t., (2.9) 

t < tc, 

where t = log[(-j92 - g^)/yljj, tc = -2, tp = 0.5, Aged = 592(MeV), and A^^ and Nf are 
the numbers of colors and flavors respectively. Here, we use N^ = 3 and Nf = 2. Bn{x) 
depends on the regularization parameter tp. For parameters in running coupling, we use the 
parameters in Ref. 8) (see §3.2). 

§3. Results of numerical calculation 

3.1. Functional gauge parameter 

The WTI is satisfied by taking the Landau gauge at zero temperature, as shown in §2.11 
By contrast, the WTI is not satisfied using a constant gauge parameter at finite temperature. 
In fact, A'j^{x), Cn{x), and An{x) are not unity at finite temperature even if one takes the 
Landau gauge (see Figs. [3]and|l]). Hence, we assume that the gauge parameter is a function 
depending on the external momentum. 

We divide Eqs. (12 -Gap . (12 -Sap , and (]2-8bp into the gauge parameter term and no gauge 
parameter term, e.g., 

Cn{x) = l + ^X„{x) + Y^{x). (3-1) 

In this equation, we treat ^ as a function dependent on the external momentum Un, x. 
Since a gauge parameter exists in the gluon propagator, it should fundamentally depend on 
a momentum of gluon (see Eq. (l2-3p ). However, for simphcity, we use the functional gauge 
parameter that depends on only the external momentum. (Although we do not write in 
terms, C,n{x) clearly depends on temperature.) From this assumption, we configured that 
the gauge parameter must satisfy ^„(a;)X„(x) + Yn{x) = 0. By this method, we can perform 
numerical calculation satisfying the WTI. 

Note that there is a problem for Eqs. (12 -Sap and (l2-8bp . It is the fact that the temperature 



dependences of C„(x) and A„(x) are different. Owing to this property, it is difficult to make 
Cn{x) = 1 and An{x) = 1 at the same time. However, Cn{x) and An{x) have near values 
at around the critical temperature. Hence, in this important region, we can simultaneously 
make C„(x) ~ 1 and An{x) ~ 1. We operate C„(x) and A„(x) so that those may approach 
1 as much as possible. 



3.2. Numerical calculation 

Although an effective potential is needed to study the phase transition, B„{x) = directry 
denotes a critical point at finite temperature.'^ Thus, we focus our attention on the critical 
point at Bn{x) = 0. 

Note that our numerical calculation has an error due to the number of flavors in the 
running coupling. In Ref. 8), the running coupling with Nj = 3 was used. However, we 
took Nf = 2. This difference in Nf in the running coupling with the same Aged results in 
the difference of about 2 (MeV) for the pion decay constant. However, in this paper, we 
are designed to search for the existence of a solution satisfying the WTI and the shift of 
the critical temperature. Thus, this difference in the pion decay constant is not important 
for that purpose of the study. (In addition, this difference might be within the range of 
numerical error.) The critical temperature has the error range of about ±3 (MeV) at least 
by ignoring this. 

We calculated the following four cases: 

(I): coupled equations, A'^{x) and Bn{x) with the Landau gauge, 

(II): coupled equations, C„(x),y4„(x) and Bn{x) with the Landau gauge, 

(la): coupled equations, A'^{x) and Bn{x) with the functional gauge parameter ^^(x), and 

(Ila): coupled equations, Cn{x), An{x) and Bn{x) with the functional gauge parameter ^„(a;). 

To solve the SDE, we employed the iteration method. It starts as a constant or a trial 
function, and repeats until a value is converged. First, we tried the trial function like zero 
temperature (see, e.g., Refs. 17) and 19)). The trial function has the zero temperature form 
for X at n = 0, —1; others are small constant values. However, there is no difference between 
a constant and the trial function except the convergence. Thus, we used a constant. 

The range of summation m = — 10~9is sufficiently large as a truncation point for 
summation in the case of A'^{x) = 1. In contrast, in (I) and (II), since the range of —10 ~ 9 
is insufficient (especially in the calculation of the critical temperature), the summation is 
necessary to take above m = —40 ~ 39. 

Since C„(x) and v4„(x) simultaneously do not make Cn{x),An{x) ^ 1 below the critical 
temperature, (Ha) is used only for determining the critical temperature. Then, we fix the 
lowest value of C„(x) as 0.94 and the highest value as 1, because it is possible to limit the 
difference from 1 to about 0.05 at around the critical temperature. 



If one fixes C„(x) = 1 exactly at all regions, A„(x) is further away from 1, compared 
with the case of C„(x) = 0.94 ~ 1. Thus, we do not fix C„(x) = 1 exactly even around the 
critical temperature. We fix A'^{x) = 1 for (la) and limit C„(a;) = 0.94 ~ 1 for (Ila). 

On the other hand, we assume that C„(x) = An{x) = 1 in (Ila) has a problem, because 
C{p) 7^ A(p) {po is continuous) in real time generates plasminos for fermion.'^SI' From this 
viewpoint, C„(x) ^ An{x) ^ 1 might be valid as an approximation (see also §4). 

The convergence for -Bo(O), Co(0) and Aq{0) in (II) is shown in Figs. [Hand [21 The conver- 
gence is worse around the critical temperature. Other cases also have the same convergence. 

3.3. Results of (I) and (II) 

The typical x and n dependences of v4^(x) and C„(x), An{x) are shown in Figs. [3] and 
[H (Figure [3] shows only the case of A^(0). The behavior for C„(0), A„(0) is much the 
same.) The temperature dependences of Ao(0) and Co(0), Aq{0) are also shown in Fig. [3 
Those results actually show A'^{x) ^ 1 and C„(x) ^ 1, v4„(x) ^ 1. In particular, Aq(0) and 
Co(0), Ao{0) shift from 1 mostly in the vicinity of the critical temperature. (The critical 
temperature is shown in Fig. [T2J) Thus, A'^{x) and C„(x), A„(x) strongly contribute to the 
critical temperature. The effect of A'^{x) or C„(x), An{x) on the critical temperature with 
A'^{x) = 1 is estimated at about 30 (MeV). 

3.4. Results of (la) and (Ha) 

In (la), A'j^{x) is unity at all regions. In (Ila), the behaviors of Cn{x) and An{x) with 
the functional gauge parameter are shown in Figs. [H] and [3 The dependence of x for n 7^ 
is also a similar form. Owing to forcibly making the situation C„ ~ 1, n and x dependences 
of An{x) change from the case of (II). Nevertheless, since n and x dependences of Bn{x) are 
unchanged in the result of numerical calculation {Bn{x) is shown, e.g., Refs. 8) and 9)), we 
expect that the properties of Bn{x) are not lost. 
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Fig. 4. X dependences of Aq (a;) and Co (a;) , Ao{x) 
at r = 0.12 (GeV). 
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Fig. 5. Temperature dependences of Aq{0) and Co(0), ^o(O). Those have a peak around the 
critical temperature. 



Cn{x) and An{x) have values that are closer to 1 than the result with the Landau gauge, 
e.g., Co(0) = 1.264 and Ao{0) = 1.368 in (II), Co(0) = 0.943 and Ao(0) = 1.139 in (Ila) at 
T = 0.12 (GeV). Moreover, C„(x) and An{x) approach 1 simultaneously above the critical 
temperature. (Co(0) and Aq{0) are shown Fig. [SI Other cases also have this feature). How- 
ever, as mentioned in §3.1, this method cannot achieve C„ ~ 1 and An — 1 simultaneously 
at low temperature. 

The functional gauge parameters ^^(x) and C,n{x) are shown in Figs. [9]andfT0l ^o-i(^) 
(^0,-1(3^)) has the most different value from the Landau gauge ^ = 0. 

^^(x) and ^„(x) have similar values at around the critical temperature, because A'^{x) 
resembles Cn{x) to some degree (see Figs. [Hand [5]). However, there is a different behavior 
above x = 1 (GeV) (Fig. [9]). This difference is understood from Eq. (13 -ip and Fig. |H In Eq. 
(l3Tp . if Cn{x) with the Landau gauge is below 1 (F„(x) is negative) and X„(x) is positive, 
^n{x) is positive. Similarly, if C„(a;) with the Landau gauge is above 1 and Xn{x) is positive. 
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Fig. 8. Temperature dependence of Co(0), Aq{0) with the functional gauge parameter. 

C,n{x) is negative etc. Therefore, the difference above x = 1 (GeV) in Fig. [H] results from a 
sign of Xn{x), Cn{x) — 1 and A'^{x) — 1. 

For large \n\, Cni^) ^^^ ^n{x) have large values. For small \n\, they have large values at 
large x. Since A'^{x) and Cn{x), An{x) with the Landau gauge are about unity at that re- 
gion, a gauge parameter does not contribute to A'^{x) and C„(x), An{x) if a gauge parameter 
is not large. Hence, the large values of ^^(x) and ^n(a^) at that region are not meaningful. 
Moreover, since A'^{x) and C„(a;), An{x) at that region are about unity, there is no problem 
in the use of the Landau gauge at that region. For this reason, it is possible to choose a 
gauge parameter as a step function depending on the external momentum. For example: 



^n{x) = a{T)Sn,x X 



A 
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1 X < Aged 

X> Aged 

where a{T) is a function that has appropriate values for each temperature. 
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Fig. 9. X dependence of ^^{x) at T = 0.15 
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Fig. 10. n dependence of ^n(O) at T 
(GeV) 



0.15 



0.05 


X 


















n (lla) " 


0.10 


D 


X 
















X (la) ■ 


0.15 




D 
















- 


0.20 






X 














- 


0.25 






D 


X 












- 


0.30 




















■ 


0.35 










X 






x 

D 


g 


s 


0.40 
0.45 










□ 

X 

nXn 


X 

D 


x" 






- 


0.50 










D 










- 



0.10 



0.15 



0.20 



0.25 
T(GeV) 



0.30 



0.35 



0.40 



Fig. 11. Temperature dependence of ^o(O), Co(0)- Co(0) ^'^^ Co(0) have the downward peak corre- 
sponding to the largest values of Co(0) and ^o(O) (^^^ ^^?>- I5D- 

3.5. Critical temperature 

The temperature dependence of i?o(0) is shown in Fig. [T21 (I) and (II) have similar 
temperature dependence and critical temperature. Their critical temperature is 141 (MeV). 
As the result, one finds that the difference in (I) and (II) hardly affects S„(x). 

In Ref. 21), the case of A'^{x) = 1 with Landau gauge and (II) was calculated with a 
different running coupling constant (Higashijima-Miransky type). This paper also showed 
that (II) has a lower critical temperature than in the case of A'^{x) = 1 with Landau gauge. 
The difference in the running coupling constant is the method of infrared cutoff and, perhaps, 
the maximum value. These differences do not strongly affect the behavior of i?„(x). Thus, 
the temperature dependence of Bn{x) for these two cases is similar to this paper. 

The critical temperature of (la) is 160 (MeV), (Ila) is 161 (MeV). Those also have almost 
the same value. The critical values of (la) and (Ila) are about 10 (MeV) lower than the result 
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Fig. 12. Temperature dependence of Bo{0) in the cases of (I), (II), (la), (Ila), and the Landau 
gauge with A'^{x) = 1. 

of the Landau gauge with A'^{x) = 1. 

On the other hand, since the choice of C„(x) contributes to Bn{x), the critical temper- 
ature is also affected to some degree. If we fix Cn{x) = 1, the critical temperature is 159 
(MeV). The result in Fig. [12] is the case in which Cn{x) and v4„(x) are simultaneously near 
1 within 10~^ order. 

Consequently, C„(x) = An{x) in the general form of the exact quark propagator is a 
reasonable approximation. Moreover, if (Ila) is more correct from the viewpoint of gauge 
invariance, the SDE with the Landau gauge and A'^{x) = 1 is the simplest, reasonable 
approximation in the framework of the ladder approximation. 

§4. Summary and discussion 

The improved ladder approximation SDE at finite temperature has a gauge choice prob- 
lem. For this problem, we calculated the SDE with the gauge parameter depending on an 
external momentum. By this method, the WTI is satisfied. Then, we used two cases for the 
exact quark propagator form, that is, C„(x) = A„(x) and C„(x) 7^ An{x). 

The result with the Landau gauge for the critical temperature shows that C„(x) = An{x) 
is a reasonable approximation for the general form C„(x) 7^ A„(x). Thus, when one studies 
the critical temperature using the SDE, C„(a;) = v4„(x) is the valid method. 

To solve the SDE with the functional gauge parameter, we found that the functional 
gauge parameter has essentially large values for, e.g., n = — 2 ~ 1 and small x(< Aged)- 
Thus, we expect that it is valid to choose the functional gauge parameter as a step function 
depending on the external momentum. 
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In the case with the functional gauge parameter, the critical temperature is nearly the 
Landau gauge with A'^{x) = 1. For this reason, we found that the Landau gauge with 
A'^{x) = 1 is a reasonable approximation from the viewpoint of the WTL 

Finally, we point out the uncertain part with respect to thermal effects in our calculational 
procedure. C„(x) and An{x) include thermal effects. In particular, for real time, C{p) ^ A{p) 
(po is continuous) has a physical meaning. This generates characteristic collective modes, 
plasminos (for fermion). Thus, rigorous C{p) = A{p) = 1 might have a problem. Therefore, 
for real time, it is difficult to satisfy the ladder approximation WTI fully. In view of this, 
Cn{x) ^ An{x) by numerical method in our calculation might be valid as an approximation. 

On the other hand, we assume that this problem might not be critical in our calculational 
procedure with imaginary time, because the analytic continuation is necessary to change from 
imaginary time to real time, after performing the summation. Moreover, the relation between 
imaginary time and real time is not simple. For example, as shown in Ref. 7), there is an 
extra term in the SDE with real time. Hence, our calculational procedure with imaginary 
time should not provide correct results in real time. If we can make C„(x) = An{x) = 1 
from our calculational procedure, it should become C{p) = A{p) = 1 after the analytic 
continuation. However, we do not know whether this should be meaningful in real time. At 
least, we must study the existence of a solution corresponding to C{p) = A{p) = 1 after the 
analytic continuation. In addition, since the functional gauge parameter that we used here 
depends on an external momentum, we must modify the functional gauge parameter in some 
way in real time. Therefore, we hypothesize that our method does not affect the existence 
of plasminos directly. 

To understand the exact details, we should study the relation between real time and 
imaginary time by the same method used in Ref. 7). 

Appendix A 

L, I and H in the SDE 



We show explicit expressions of L, I, and H in Eqs. (12 -Gap . (12 -Sap , and (]2-8bp . 
a+ = {po - qaf - {x + yf , a_ = (po - qof - {x - yf. 



A'ix) 



x^ + y^-pl- ql a+ 

Li = log — + 2xy, 

2 a_ 



(Po -qo-x^ + y ) fJ: L 
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Cn{x) 



AJx) 



h = 2go log — , /2 = go ( - 2(po - qofy 



h = -(Po -go) 



log — 
a_ 



i 

a+ 


-^) 


log 


-x' 


-.M 


.a_i_ a_/J 



Hi = (po - qo)qo 



[X' -y' + (po - go)') {- ) + log — 



H2 = -Axy + (x' + y'- {po - go)') lo, 



H. 



X +1/ ^ jlog- 



{x^ -y^Y - (po-go)S / 1 
2 y Va+ 
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